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MULTIPLE POTENTIAL WELL STRUCTURE IN INERTIAL ELECTROSTATIC 
CONFINEMENT DEVICES 


Ryan M. Meyer 

Dr. Sudarshan Loyalka, Thesis Supervisor 
Dr. Mark Prelas, Thesis Supervisor 

ABSTRACT 

Inertial Electrostatic Confinement (IEC) devices are of interest as neutron 
generators for many applications. An essential part of these devices is the formation of a 
multiple potential well structure within the devices. In this paper, previous analyses of 
the formation of these wells are reviewed and extended. Three types of IEC systems are 
classified and analyzed according to the arrangement of electrodes and the species within 
the system. These systems are the uni-polar cathode-anode (UCA) system, the bi-polar 
cathode-anode (BCA) system, and the bi-polar anode-cathode (BAC) system. Results of 
extensive parametric studies are reported through an efficient solution of the Poisson’s 
equation. These results delineate the conditions most conducive for double potential well 
formation in different systems and may aid in the future design of IEC systems. 
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INTRODUCTION 


Inertial Electrostatic Confinement (IEC) devices have many potential 
applications, including their use as neutron generators. A few of these applications 
include: luggage inspection, oil well logging, medical isotope production, detection of 
explosives, archaeometry, forensics, breeding advanced fuels (He3), transmuting long- 
lived radioactive isotopes from fission reactor waste, a neutron source to drive a sub- 
critical fission reactor, generating power for space craft, and generating power for 
terrestrial needs [l]-[9]. 

Medical isotopes are typically produced on site at a clinic with the use of medical 
isotope generators, or accelerators. Isotopes used for imaging must be produced on site at 
the clinic due to the short half-life of typical isotopes employed for medical imaging. 
The advantages of an IEC device for producing medical isotopes are that it is small, 
compact, and inexpensive. Such a portable neutron generator can be transported to some 
of the most remote places on earth in order to treat patients in a relatively inexpensive 
manner. Another advantage of the IEC is that the radioactivity hazard can effectively be 
“unplugged”, when not in use, thus, special handling, and shielding requirements are 
limited [10]. Of course, many of these characteristics are desirable features for a number 
of neutron generator applications. Experiments performed at the University of 
Wisconsin-Madison have generated 4-8 Bq of N13. Protons of 14.7 MeV, generated by 
D-He3 reactions, irradiate water flowing through the IEC reactor to produce N13. The 
isotope, N13, is later separated from the water [7], [8]. 
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In light of these potential applications for the IEC devices, many experimental 
and theoretical studies have been conducted regarding the potential distribution within 
IEC reactors. Experiments performed with IEC devices indicate the existence of a 
multiple potential well structure, within the inner electrode. The general method of 
theoretically obtaining the potential distribution within a spherical or cylindrical IEC 
device has been to solve Poisson’s equation while simultaneously applying the conditions 
of current continuity and conservation of energy. These calculations, however, have not 
consistently supported the formation of double well structures, indicated by experiments 
[1 1]-[15] . It is, therefore, of interest to review the previous theoretical work, and to study 
the well formation in one consistent framework of computations. 

Although Langmuir and Blodgett [16] did not study the IEC devices explicitly, 
their work still provided the basic solutions for potential distribution in the systems of 
interest. These authors solved Poisson’s equation within the inner electrode of a system 
of concentric spherical electrodes, in a vacuum: 


ii l2 ^ = i r l'2. 

X 2 dX dX X 2 


( 1 ) 


With the boundary conditions for space charge limited current [17] (i.e. — — = 0 and 

dr 

V =0 at the emitting electrode). In Eq. (1), Y is the normalized potential and X is the 
normalized radius, given by: 


f= m 

V a 


( 2 ) 




r 

~R a 


( 3 ) 


While K e is defined by Eq. (4), where I e is the electron current, m e is the electron mass, 
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( 4 ) 


e is the magnitude of the electron charge, while V a and R a are the anode voltage and 
radius, respectively. The typical solution generated by Langmuir’s and Blodgett’s 
equation, Eq. (1), is provided in Fig. 1 for the region within an inner anode. 
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Figure 1: Langmuir’s and Blodgett’s solution to Poisson’s 
Equation for concentric spherical electrodes [16]. 

Hirsch [11] considered a system involving ion injection into a grounded vacuum 
chamber and electron emission from the inner surface of the chamber wall. For his 
model, Hirsch constructed the equation, for the region within an inner cathode, 


J__ d _ x 2 — 

X 2 dX dX 


|l [r-^-^a-n-' 2 ]; 


(5) 


along with the boundary conditions for space charge limited current [17]: 


Y ( 1 ) = 0 , 


dX 


( 6 ) 


Here , and K t are defined by 


4 = 




1 

47T£ 0 


(7) 

( 8 ) 
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and the normalized potential is redefined in Eq (9), as: 


Y -~V(r) 

V c 


( 9 ) 


In Eqs. (7-9), 7 ; is the ion current, m ( is the ion mass and - V c is the cathode potential. 
Hirsch found that the solution to the boundary value problem, defined by Eqs. (5) and (6), 
consists of a double potential well structure within the inner cathode, as shown in Fig. 2. 

Potential vs Radius 
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Figure 2: Hirsch’s solution to Poisson’s Equation for ions injected from a spherical vacuum 
chamber, with electrons emitted from the inner surface of the vacuum chamber [11], 


Hirsch’s model neglects spreads in the total and angular energies of the particles. 
However, Dolan, [15], [18] Lavrent’yev, [19] and Swanson, [13]-[14] created more 
realistic models that allow for spreads in the total and angular energies of particles. 
Dolan and Swanson solved Eqs. (10), and (1 1) 

_L J_ x 2 — = f(X,Y) (10) 

X 2 dX dX Je 

J_ ^_x 2 — = f e (X,Y)- gi (X,Y) (11) 

X 2 dX dX e 1 

where f e (X,Y), and g t (X,Y) are given as 
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fe(X,Y) = 


v a 


4 7T£ 0 { 

X 2 (Y-Y e>E 


( 12 ) 


JX* \dr , 4 


l e (Ye,E^) 

X 2 [Y-Y eE -(Y e JX 2 )f 


8i (X,Y) = - 


4 7ie 0 { m A 

1 X 2 (Y iE -Y) 


(13) 


\ dY TJy 

r “ l X 1 V, e -Y-(Y l JW 


for a system with electron emission from an outer grounded cathode and for a system 
which included ions confined within the inner anode as well as electrons emitted from the 
outer cathode. The normalized total and angular electron energies are represented by Y e E 


and Y e ^ , while the norm a li z ed total and angular ion energy is represented by Y t E and 
Y i4 . Equation (10) produces broad single well solutions, while Eq. (11) reveals the 
formation of double wells under certain conditions. 

Swanson’s results [13] indicated that double well solutions do form under limited 
circumstances, however, questions arose pertaining to the uniqueness of the double well 
solutions. Dynamic IEC analysis by Hockney [20] and Ohnishi [21] indicated that 
double potential well structures are transient phenomena, thus double well solutions 
computed in the steady state may not be valid. In addition, simulations by Matsuura [22] 
revealed that a multiple potential well structure is not necessary to produce multiple 
radial peaks in the neutron production rate. 

In the future, for clarity in system designation, IEC systems will be denoted as 
shown in Table I: 
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Table I: This table identifies different types of inertial electrostatic devices that will be analyzed in this work. 


Denotation Electrons Ions Schematic 



Figure 3: Schematic diagram of a UCA 



Figure 4: Schematic diagram of a BCA 



Thus, the studies performed by Langmuir and Blodgett were for a UCA system, while 
Hirsch’s model is equivalent to the BAC system. In addition, Eq. (10) represents a UCA 
system and Eq. (1 1) is for a BCA system. 

Our purpose here is to obtain and critically examine potential well solutions 
within the inner electrodes of the systems in Table I. We also have explored further 
results for the UCA system and show that we get results similar to those of Dolan and 
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Swanson [12] -[15]. The method of analysis for the bi-polar systems will involve finding 
a double well solution and performing systematic variation of certain parameters to 
determine which parameters affect double well formations the most. We should note that 
more detailed IEC modeling is performed by simultaneously solving Poisson’s equation 
with equations describing the creation and removal of particles from the system, while 
the present paper is directed at careful review and analysis of the problem where specific 
forms of particle distributions are assumed. This work thus serves as a benchmark for 
more detailed considerations. 
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II. EQUATIONS AND METHODS OF SOLUTION 


For the systems in Table I, the basic problem can be described by Eq. (10) for the 
UCA system and by Eqs. (11) and (14) for the BCA and BAC systems. 

J__d_ x 2 — = f.(X,Y)-g e (X,Y) (14) 

X 2 dX dX ' 

InEq. (14), f.(X,Y ) , and g e (X,Y) are, respectively; 


fi(X,Y) = 




Y X\Y-Y UE ) j ( y y , 

\dY iE \dY — ' UE ' 

o \ X\Y-Y. E -(Y.JX 2 )f 


(15) 


g e (X,Y) = - 


4 K£ c 

X 2 (Y etE -Y) 


\dY eE jdY j — - — 

j e ’ J e4 X 2 \ Y „-Y-(Y Jx 


(16) 


^ X 2 [Y e E - Y - (Y e ^ jx 2 )] 1/2 
If we assume guassian distributions for current in Eqs. (12) and (13), then for the systems 
described by Eqs. (10) and (11), the electron current is distributed over the normalized 
electron total and angular energies, Y e E and Y e ,, as shown in Fig. 6. In addition, the ion 
current is distributed over T E and Y i(j) , as shown in Fig. 7, where Y i E and T ^ are the 


normalized ion total and angular energies, respectively. 





0 1 

Figure 7: Gaussian ion current distribution over the Y i E , plane. 

The electron current may be represented with Eq. (17), 

= (17) 

with cr eE and o e 0 representing the spreads in total and angular electron energy, and 
where the constant A e is computed by recognizing that Eq. (18) represents the circulating 
electron current. 
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Ice=V«\dY e , E jdY e / e (Y eE ,Y e ^) 


(18) 


One thus [13], [14] has for A e : 

t - e 4/0 - [<j e , l(a e ^ - o e E )\ [e~ iy ^ ] - 

Evaluating A,, is not as straightforward as the evaluation of A e because the electron is 
allowed to have Q<Y e(/) <1-F c£ and 0 < Y e f < 1 . The ion has energies in the range 
0 <Y i(/) <{y iE - T min ) and y mjn <Y iE <\ . These limits depend on the minimum norm a li z ed 
potential within the anode, Y mm , which is information that we seek from the solution of 
Eq. (11). The ion current is represented with Eq. (20), and the circulating ion current is 
defined in Eq. (21). 


I,v,„ = 

(20) 

l 

i cl | dY iE 

(21) 


Observation of Fig. 7 reveals that for small values of energy spreads in the total 
ion energy and angular ion energy, <7 i E and <7,.^, integration of Eq. (21) may be 
approximated with T min =0 for moderate values of Y mm (i.e. Y mm < 0.85), which allows 
us to compute A, from Eq. (22). 

I- 
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If Uj E =0.05 and cr ^ = 0.00 1 , then values of A i , for values of Y mm from 0 to 1 are 
provided in Table II. 


Table II: A, versus Y mm for a i>E = 0.05 
and 0^ = 0.001. 



A/(Ii V a ' 2 ) 

0 

20000 

0.1 

20000 

0.2 

20000 

0.3 

20000 

0.4 

20000.1 

0.5 

20000.9 

0.6 

20006.8 

0.7 

20050.7 

0.8 

20380.9 

0.85 

21070.4 

0.9 

23204.5 

0.95 

32019.8 


The error in using the approximation T mm =0 when T min =0.85 is 5.1% and the error in 
using Y mm = Owhen T min < 0.80 is < 1.9% . Using Eqs. (17) - (22), we can introduce Eqs. 
(23) and (24), 


f e (X,Y) = 
K 


a f x2{¥ r ¥ ^ 

x^[ dYeE [ dY ^vY-y e , E -Klx 2 f 


8i(X,Y) = 
K. 


a r x( Y r) 

-^f dY iE \dY u — , ™ 

X 2 { ‘’ E l ^[Y, e -Y-{ Y jX 2 )! 1 ' 


(23) 


(24) 


with K e and K t defined in Eqs. (25) and (26), where /; represents the transparency of the 
anode [15], [18], 
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( 25 ) 


K e = Mi-t 2 )]- {[ i/fe ^.^J-k/NF}- 

(avr) 

K t =k/(l -7 2 )]- {[l/(4^ o a iE o;.J]-k/(2,)] 1/2 }. 

{i-^ + k,/k,-^j-k^^]r- ( 26 ) 

(i./vr) 


Eqs. (23) and (24) may then be plugged back into Eqs. (10) and (11), respectively, to 
solve for the potential distribution within the anode of a UCA system and a BCA system. 

In the BAC system, represented by Eq. (14), ions and electrons reverse the roles 
they had in the BCA system of Eq. (11). The distribution of ion current for Eq. (14) is 
now represented with Fig. 6, and the distribution of electron current is now illustrated in 
Fig. 7. The ion current is represented with Eq. (27), 

= A- (27) 

and the circulating ion current is now defined with Eq. (28). 


hi=Vc\dY LE pVA*’ 1 *) 


(28) 


A is surmised from Eq. (28) and defined in Eq. (29). 




(29) 


Electron current is represented in Eq. (30), and the circulating electron current is now 
defined by Eq. (31). 

I e (Y eE , Y e<p ) = A e AK (30) 
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Analogous to A , for the BCA system, a difficulty occurs in determining the value of 


A e , for the BAC system, because ranges for the electron total and angular energies are 
Q<Y e(j) <iy eh: - y min ) and Y nin <Y eE < 1. As in the case of the BCA system, this BAC 
system requires knowledge of the minimum normalized potential, which is information 
that we seek from the solution of Eq. (14). Thus, T min is once again approximated with 
Y mm = 0 to obtain the value of A e in Eq. (32).: 

A,=L/(v/^ CT ,. s )j- 

(i-^ ! + k,/k, -Ol 

This assumption leads to errors in the calculations of A e of -40% for large values of 
Y mm (-0.85). An error of 40% will produce a noticeable difference in the solution, but, it 
will not produce a change in the trends observed. In addition, since the error decreases 
exponentially as the value of y min decreases, error in the calculation of A e for the BAC 
system, does not produce a noticeable change in the solution for small values of T mm . A 
table with values of A e calculated using various values of Y mm is given in Table III. 



Table III: A e versus for a e<E = 
0.15 and a e ,<|,= 0.001. 


Y^n 

Ae/(Ie V c ' 2 ) 

0 

6675.22 

0.1 

6683.34 

0.2 

6699.23 

0.3 

6730.38 

0.4 

6791.9 

0.5 

6915.01 

0.6 

7168.07 

0.7 

7718.23 

0.8 

9074.82 

0.85 

10587.9 

0.9 

13798.7 

0.95 

23924.0 


In parallel with Eqs. (23) and (24), for BCA systems, we obtain Eqs. (33) and (34) to be 
plugged into Eq. (14). 


The values of K t 
and (36), where rj 


fi(X,Y) = 


Ay,,eKe)Ay^K 


x 2 ]<U, \ B jdY„ 

g,(X,Y) = 


far.,- 


AYe,El<Te,E)AYe,t/°e.<,) 


(33) 


(34) 


*lY e 'E-Y-(Y e JX 2 } ] 1/2 
and K e for Eqs. (33) and (34), respectively, are defined in Eqs. (35) 
is the transparency of the cathode. 

K i =Mi- 7 2 )]-{[i/( 4 ^ J- k/N] 1/2 }- 

fi-^)_u /U 


(35) 


(1,/vr) 
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( 36 ) 


K e =Ml-7 2 )]- {[ V( 4 Ke „ h/{2e)] l/2 }- 



Equations (10), (1 1), and (14) are solved numerically. The double integrals on the 
right hand sides of Eqs. (23), (24), (33) and (34) have been previously evaluated by 
expanding the exponential into a Taylor Series. After Taylor series expansion, the 
integral of the resulting argument may be found using a standard table of integrals [13]. 
However, we have found that it is possible to express the double integrals on the right 
hand sides of Eqs. (23), (24), (33), (34) in terms of error functions for which fast and 
accurate numerical routines are readily available. Thus, for g i (Z,T)and g e (X,Y), we 
have: 

8i(X,Y) = {kJxYItt 112 a. E ajj(a i0 ~(7 , e X 2 \ 

k; 2 . . Erfux (1 _ y) 172 /^; 2 ] - | (37) 

g e (X,Y) = (. K e /X ) • \W ' 2 a e h: aJl(<J e j-<r e , E X 2 )]- 

k . . Erfi[X (1 - Y ) U2 /<j' e [ 2 ]- 1 (38) 

I o-] 12 e [{Y -' )/(Tj • X • Erfi[{\ - Y) y2 /al l2 ]\ 

where Erfi[z] is the imaginary error function, which is defined in Eq. (39), where Erf[z\ 
is the error function. 

Erfi[z\ = ^^ (39) 

i 

Similarly for f e (X,Y) and f t (X,Y) we find: 
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( 40 ) 


f e (X,Y) ={K e lx)\n m <r eE aJl((T e j -<J e , E X 2 )]- 
|<; 2 ■e^^-Erfi\XY l,2 /o^]- 
I <1 ■X-Erfi[Y l,2 lal l2 \ 


ft(X,Y) = {KJX)- [k 2 cx., aJl(a i 0 -a hE X 2 \ 

2 . k rfa i . Erfi [X y V 2 /(t} 12 ] - 

o?je^ ) -X>Eifi[Y U2 /ot£] 


(41) 


These results dramatically reduce computational labor and increase computational speed. 
In addition, they increase the accuracy and reliability of solutions obtained. 

Equations (37), (38), (40), and (41) can be inserted into Eqs. (10), (11) and (14). 
The latter equations are then solved numerically by assuming the boundary condition in 


Eq. (42), 


dY( 0) _ Q 
dX 


(42) 


and employing the shooting method to determine the value of 7(0) that satisfies the 
boundary condition of Eq. (43). 

7(1) = 1 (43) 


The models just described, for each of the systems in Table I, are summarized in Table 
IV. 


Table IV: A summary of the systems and corresponding equations. 


System 



Model 

Right Hand Side of Model 

UCA 

1 d X 2 

■ dY 

= f e (X,Y) 

Use f e {X,Y ) inEq. (40), with K e inEq. (25). 

V 2 dX ' 

dX 

BCA 


■ dY 

= f e (X,Y)- gi (X,Y ) 

•Use f e {X,Y ) inEq. (40), with K e inEq. (25). 

V 2 dX 

dX 

• Use g t (X,Y) inEq. (37) with K t inEq. (26). 

BAC 


■ dY 

= MX,Y)-g e (X,Y ) 

•Use f\{X,Y) inEq. (41), with K i inEq. (35). 

V 2 dX 

dX 

• Use g e {X,Y) inEq. (38) with K e inEq. (36). 
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In order to study double well solutions for Eqs. (11) and (14), it is convenient to 
introduce some measures. A central well will provide better confinement if the peak of 
the well is large and the well is narrow in width. The width of a central double well is 
depicted in Figure 8, 



Figure 8: Depiction of a typical double well solution. 


and is defined as the Full Width at Half Maximum (FWHM). The value of half- 
maximum is the average of the central well peak potential and the minimum potential 
within the well. The “Double Well Depth” (DWD) is defined by Tzonev [22] in Eq. (44), 
and is depicted in Figure 8. 

Y -Y 

DWD = — (44) 

1 - Y mm 
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III. RESULTS 


Solutions for systems considered in Table IV will be provided in this section. A 
parameter for each system will be varied, while the rest of the parameters are held 
constant. Thus, several plots are generated for each of the models in Table IV, revealing 
trends in IEC behavior. The FWHM and DWD of solutions in a particular plot are 
provided in a corresponding table. We begin by discussing the results for the UCA 
system, and then provide the results for the BCA and BAC models. 

A. UCA System: 

The solutions to Eq. (10) for different values of o e E , o e ^ , I e , and V a are shown 
in Figs. 9-12. Figure 9 shows several solutions to the UCA model for values of <T e ^ over 
the range of 0.0005<<J C ^<0.8 . Solutions, for multiple values of cr eE , where 
0.01<<T C E <0.8 are shown in Fig. 10. Figure 1 1 contains solutions for UCA model as the 
voltage is varied from 50<V a <2000 V, while Fig. 12 contains solutions for several 
values of current, with 25 <I e <375 mA. 

Figure 9 shows that the solution is highly dependent on the angular energy spread, 
and, as the angular energy spread of electrons decreases, the solution of Eq. (10) 
approaches the solution of Langmuir and Blodgett [16] in Fig. 1. Variation in the total 
energy spread of electrons has very little effect on the potential well shape as indicated in 
Fig. 10. Figures 11 and 12 show that potential well depth varies proportionally to the 
1 

perveance, — fjj . 
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Potential vs Radius 
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Figure 9: Potential vs. Radius within the anode of the 
UCA model as a function of angular energy spread. 



Figure 10: Potential vs. Radius within the anode of 
the UCA model as a function of total energy spread. 


Potential vs Radius 



Figure 1 1 : Potential vs. Radius within the anode of 
the UCA model as a function of applied voltage. 



of the UCA model as a function of current. 


19 


B. BCA System: 


Each one of the parameters, <7 e ^ , cr eE , <7 i<t) , <T I E , I e , V a , and /? (ratio of the ion 

current to electron current, — ), are varied, while the rest are held constant to observe 
I e 

trends in the solution of Eq. (11). Figure 13 shows solutions for the BCA system for 
various values of <7 i</) in the range 0.001 <cr^ <0.5 , and Table V provides a list of their 
FWHM’s and DWD’s. Fig.14 contains several solutions for o iE over 0.0\<<r i E <0.8 , 
with the FWHM and DWD of each solution listed in Table VI. Electron angular energy 
is varied from 0.01<<T C ^ <0.5 and solutions for these values are given in Fig. 15, with the 
associated FWHM and DWD data provided in Table VII. Next, solutions for <j eE over 
0.05<(7 eE <0.S are shown in Fig. 16, and the FWHM’s and DWD’s for these solutions 
are given in Table VIII. The ratio of ion current to electron current, /? , is varied for 
0<J3<0.003 and solutions are given in Fig. 17, while the FWHM and DWD data for 
these solutions are given in Table IX. The applied voltage is varied for 400 <V a <1500 
volts and solutions are given in Fig. 18, with a list of the FWHM’s and DWD’s given in 
Table X. Finally, Eq. (11) is solved for several values of electron current, 
100<7, < 400 mA, and the solutions are displayed in Fig. 19, with a list of the FWHM’s 
and DWD’s provided in Table XI. 
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Potential vs Radius 
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Figure 13: Potential vs. Radius within the anode of the 
BCA model as a function of ion angular energy spreads. 


Figure 14: Potential vs. Radius within the anode of the 
BCA model as a function of ion total energy spreads. 
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Figure 15: Potential vs. Radius within the anode of the BCA 
model as a function of electron angular energy spreads. 


Figure 16: Potential vs. Radius within the anode of the 
BCA model as a function of electron total energy spreads. 



Figure 17: Potential vs. Radius within the 
anode of the BCA model as a function of (3. 
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Figure 18: Potential vs. Radius within the anode 
of the BCA model as a function of voltage. 


Figure 19: Potential vs. Radius within the anode of 
the BCA model as a function of electron current. 
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Table V: DWD and FWHM for the BCA 
model, for several values of ion angular energy 
spread, as shown in Fig. 13. 


Table VI: DWD and FWHM for the 
BCA model, for several values of ion 
total energy spread, as shown in Fig. 14. 


<Ji,* Ymin Ypeak 


0.001 

0.005 

0.008 

0.010 

0.050 

0.100 

0.500 


0.860 

0.824 

0.799 

0.783 

0.672 

0.637 

0.594 


0.935 

0.833 

0.799 

0.783 

0.672 

0.637 

0.594 


DWD FWHM 

0.538 0.313 

0.053 0.239 


. .4». Ywfc . 

0.010 0.833 0.940 

0.050 0.860 0.935 

0.150 0.800 0.855 

0.500 0.729 0.756 

0.800 0.713 0.733 


DWD FWHM 

0.641 0.271 

0.538 0.313 

0.272 0.242 

0.099 0.184 

0.072 0.170 


Table VII: DWD and FWHM for the BCA 
model, for several values of electron 
angular energy spread, as shown in Fig. 15. 



Y™ 

Ypeak 

DWD 

FWHM 

0.010 

0.122 

0.150 

0.032 

0.077 

0.050 

0.557 

0.610 

0.121 

0.127 

0.100 

0.751 

0.832 

0.326 

0.204 

0.500 

0.959 

0.984 

0.617 

0.610 

Table IX: DWD and FWHM for the BCA 

model for several values of (3, as 
Fig. 17. 

: shown in 

3 

Y™ 

Ypeak 

DWD 

FWHM 

0 

0.582 

0.582 



0.0005 

0.632 

0.632 



0.001 

0.669 

0.685 

0.050 

0.132 

0.003 

0.860 

0.935 

0.538 

0.313 


Table VIII: DWD and FWHM for the BCA 
model, for several values of electron total 
energy spread, as shown in Fig. 16. 


CJ e ,E 

Yrfn 

Ypeak 

DWD 

FWHM 

0.050 

0.860 

0.935 

0.538 

0.313 

0.100 

0.857 

0.935 

0.545 

0.314 

0.500 

0.858 

0.942 

0.591 

0.341 

0.800 

0.862 

0.945 

0.603 

0.352 

Table X: 

DWD and FWHM for the 


BCA model for several values of 


applied voltage, as shown : 

in Fig. 18. 


V a 

(volts) 

Yrfn 

Ypeak 

DWD 

FWHM 

400 

0.196 

0.595 

0.496 

0.173 

500 

0.345 

0.626 

0.429 

0.173 

600 

0.498 

0.715 

0.431 

0.185 

800 

0.730 

0.870 

0.518 

0.240 

1000 

0.860 

0.935 

0.538 

0.313 

1500 

0.936 

0.956 

0.310 

0.335 


Table XI: DWD and FWHM for the BCA 
model for several values electron current, 
as shown in Fig. 19. 


I e (mA) Yaka Ypeak DWD FWHM 

100 0.860 0.935 0.538 0.313 

200 0.540 0.742 0.439 0.191 

300 0.315 0.613 0.435 0.172 

400 0,191 0.596 0.500 0.173 
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Figures 13-19 and Tables V-XI show that a double well solution depends most 
greatly on the value of cr i4 . A double well structure exists for small values of (J i4 , and 
vanishes rapidly as <7, ^ is increased. Small values of <J e ^ tend to deepen the minimum 
depth of the potential well, but small values of a e ^ do not necessarily enhance the 
distinctness of the central double well. Results in Table VII shows that that for 
<7^ =0.01, the FWHM of the solution is smallest for solutions in Figure 15, but its DWD 
is rather poor. On the other hand, the solution for <J e ^ = 0.5 has the best DWD value for 
solutions of Figure 15, but it has a very wide FWHM. As a general observation, it 
appears to be undesirable to have o eJj) too small, possibly because the existence of more 
electrons in the center of the device may neutralize the positive ions trapped in the double 
well. 

The value of o iE has a noticeable effect on the central double well, but the 
double well solution does not vanish nearly as rapidly for increasing values of o i E as it 
does for increasing values of <r, ^ . In fact, double well solutions still exist for relatively 
large values of cr i E . The existence of a double well appears to be nearly independent of 
< 7 eE according to the solutions in Fig. 16, while Figs. 18 and 19 reveal that the 

distinctness of the double well solution is proportional to the system perveance, . 

C. BAC System: 

Each one of the parameters, o e4 , o eE , G i (j) , o i E , /, , V c , and , is varied, with 
the rest held constant to observe trends in the solution of Eq. (14) for the BAC system. 
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Figure 20 shows solutions for the B AC system for various values of <J e ^ in the range 
0.0005 <<7 e ^ <0.1 and Figure 21 contains several solutions for a eE over 0.0 1 < a eE <0.5 . 
Ion angular energy is varied from 0.06 <o i ^ <0.2 and solutions for these values are given 
in Fig. 22. Next, solutions for <7 I E over 0.01 <(7, E <0.3 are shown in Fig. 23. (3 is 
varied for 0.1 25 </?<!. 25 and solutions are given in Fig. 24, and the applied voltage is 
varied for 500 <V C <1200 V with solutions given in Fig. 25. Finally, Eq. (14) is solved 
for many values of ion current, 2 </ ; <7.2 mA, and the solutions are displayed in Fig. 26. 
FWHM and DWD data for the solutions in Figs. 20-26 are provided in Tables XII 
through XVIII. 

In Fig. 20, it is evident that the existence of a double well structure depends rather 
strongly on, <j e ^ , since the central peak in Fig. 20 vanishes rapidly as the angular energy 
of the electrons increases. Well depth increases substantially for decreasing values of 
<7,. ^ , in Fig. 22, but the DWD of the solution is also degraded. A double well solution is 
unable to materialize for <7 eE =0.01 and <7 eE =0.05. However, Fig. 21 shows that a 
double well solution appears at approximately <r eE = 0.1, and the quality of the well then 
degrades for increasing values of o eE . Fig. 23 reveals that the value of o i E has minimal 
effects on the potential well structure, while Figs. 25 and 26 reveal that the distinctness of 
a double well varies with proportion to the perveance, —jjy ■ 
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Potential vs Radius 




Figure 20: Potential vs. Radius within the cathode of the BAC 
model as a function of spread in angular electron energy. 


Potential vs Radius 


Figure 21: Potential vs. Radius within the cathode of the 
BAC model as a function of spread in total electron energy. 


Potential vs Radius 



Figure 22: Potential vs. Radius within the cathode of the 
BAC model as a function of spread in angular ion energy. 


Figure 23: Potential vs. Radius within the cathode of the 
BAC model as a function of spread in total ion energy. 
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Figure 24: Potential vs. Radius within the 
cathode of the BAC model as a function of (3. 



Figure 26: Potential vs. Radius within the cathode 
of the BAC system as a function of ion current. 


Potentialvs Radius 



Figure 25: Potential vs. Radius within the cathode of 
the BAC model as a function of applied voltage. 
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Table XII: DWD and FWHM for the BAC 
model for several values of electron 
angular energy spread, as shown in Fig. 20. 


Table XIII: DWD and FWHM for the BAC 
model for several values of electron total 
energy spread, as shown in Fig. 21. 


Pe.t Ymin Ypeak 

0.0005 0.508 0.508 

0.001 0.760 0.845 

0.005 0.732 0.743 

0.010 0.699 0.699 

0.050 0.598 0.598 

0.100 0.564 0.564 


DWD FWHM 


0.353 0.238 

0.042 0.213 


OeJ 

0.01 

0.05 

0.10 

0.15 

0.30 

0.50 


Ymin Y pe a k DWD FWHM 


0.508 

0.508 

0.760 

0.741 

0.700 

0.674 


0.508 

0.508 

0.845 

0.812 

0.749 

0.710 


0.353 

0.275 

0.163 

0.112 


0.238 

0.224 

0.197 

0.181 


Table XIV: DWD and FWHM for the 
BAC model for several values of ion 
angular energy spread, as shown in Fig. 22. 


Pi.* 

Ymin 

Ypeak 

DWD 

FWHM 

0.06 

0.491 

0.555 

0.125 

0.132 

0.10 

0.650 

0.730 

0.229 

0.179 

0.15 

0.760 

0.845 

0.353 

0.238 

0.20 

0.820 

0.898 

0.434 

0.290 

Table XVI: DWD and FWHM for the 
BAC model for several values of (3, as 


shown in Fig. 24. 


P 

Ymin 

Ypeak 

DWD 

FWHM 

1.250 

0.554 

0.557 

0.006 

0.083 

1.000 

0.565 

0.569 

0.011 

0.096 

0.500 

0.613 

0.633 

0.052 

0.137 

0.200 

0.741 

0.812 

0.275 

0.224 

0.125 

0.836 

0.915 

0.479 

0.313 


Table XV: DWD and FWHM for the 
BAC model for several values of ion 
total energy spread, as shown in Fig. 23. 


Pi.E 

Ymin 

Ypeak 

DWD 

FWHM 

0.01 

0.741 

0.811 

0.271 

0.222 

0.05 

0.741 

0.812 

0.275 

0.224 

0.10 

0.741 

0.814 

0.282 

0.228 

0.30 

0.757 

0.832 

0.310 

0.245 

Table XVII: 

DWD and FWHM for 

the BAC model for several values of 

applied voltage, as shown in Fig. 25. 

Ve 

(volts) 

Ymin 

Ypeak 

DWD 

FWHM 

500 

0.258 

0.557 

0.404 

0.181 

750 

0.577 

0.720 

0.337 

0.203 

1000 

0.741 

0.812 

0.275 

0.224 

1200 

0.805 

0.847 

0.217 

0.230 


Table XVIII: DWD and FWHM for the 
BAC model for several values of ion 
current , as shown in Fig. 26. 

Ij (mA) Ymin Yped. DWD FWHM 

2.0 0.816 0.853 0.203 0.230 

2.8 0.741 0.812 0.275 0.224 

5.0 0.501 0.675 0.349 0.195 

1.2 0.303 0.572 0.386 0.181 
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IV. CONCLUSIONS 


Langmuir and Blodgett solved poisson’s equation for concentric spherical 
electrodes and discovered the formation of a single virtual electrode at the center 
electrode [16]. Hirsch performed similar analysis on an IEC device, but included the 
emission of electrons from the inner surface of the vacuum chamber in his theoretical 
model. Hirsch’ s experiments and analysis indicated that a multiple potential well 
structure existed at the center of the IEC device [11]. Dolan [15], [19], Swanson [13]- 
[14], and Lavrent’yev [18] developed more realistic models that accounted for spreads in 
total particle energy and angular particle energy. Following the work of Dolan and 
Swanson, we have performed computer simulations for the three systems noted in Table 
IV. We have obtained solutions in a consistent framework, and performed an extensive 
parametric study of the solutions over a e4 , a eE , <j i (/) , a iE , I e , / , V a , V c and the factor 

P- 

Information obtained, from solutions to the models in Table IV, provides insights 
in IEC performance with respect to which parameters and conditions will be most 
conducive to double well formation. The results have shown that double well formation 
requires that <J, ^ for BCA systems and o e4 for BAC systems, is relatively small. 
However, it is not necessarily desirable for o e 0 or < 7 i(Z > to be small. In that case, the 
depth of the potential well, 1 - Y mm , is increased, but the DWD, Eq. (44), is decreased. It 
could also be noted that double well existence depends weakly on <7, E for BCA systems 
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and cr eE for BAC systems. In addition, a eE for BCA systems and o iE for BAC 
systems, have very little influence on potential well structure. 

In this paper we have concentrated on models that permit extensive parametric 
studies with respect to double well formations. The models can certainly be improved 
both in terms of underlying physical phenomena and numerical techniques of solution. 
Such improvements will permit comparisons with experimental data and more realistic 
assessments. Nevertheless, the present work clarifies the past work on the problem and 
provides many useful insights. 
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